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ABSTRACT 

We use the previously-identified 15 infrared star-cluster counterparts to X- 
ray point sources in the interacting galaxies NGC 4038/4039 (the Antennae) to 
study the relationship between total cluster mass and X-ray binary number. This 
significant population of X-Ray/IR associations allows us to perform, for the first 
time, a statistical study of X-ray point sources and their environments. We define 
a quantity, rj, relating the fraction of X-ray sources per unit mass as a function 
of cluster mass in the Antennae. We compute cluster mass by fitting spectral 
evolutionary models to K s luminosity. Considering that this method depends on 
cluster age, we use four different age distributions to explore the effects of cluster 
age on the value of r\ and find it varies by less than a factor of four. We find 
a mean value of r\ for these different distributions of i] = 1.7xl0~ 8 Mq 1 with 
<j v = 1.2xl0 -8 Mq 1 . Performing a x 2 test, we demonstrate rj could exhibit a 
positive slope, but that it depends on the assumed distribution in cluster ages. 
While the estimated uncertainties in r\ are factors of a few, we believe this is the 
first estimate made of this quantity to "order of magnitude" accuracy. We also 
compare our findings to theoretical models of open and globular cluster evolution, 
incorporating the X-ray binary fraction per cluster. 
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Introduction 



The Antennae are a pair of colliding galaxies with an unusually large number of X- 
ray point sources. High resolution X-ray observation s taken with Chandra revealed 49 new 
individual sources (IFabbiano. Zezas. fc Murray! 120011). where previous observations only in- 
dicated extended filamentary structure (IFabbiano et al.lll997l ). These X-ray sources range in 
luminosity from 10 38 — 10 40 ergs s . M ost of these are thought to be X-r ay binaries (XRBs) 
with a black hole compact companion (IFabbiano. Zezas. fc Murray! l200ll ). 



In addition to numerous X-ray sources, the Antennae contain man y bright, massive 
young clusters that are evident in both optical, HST ( Whitmore et alll999 ) and infrared (IR) 
( jBrandl et al.ll2005l henceforth Paper I) images. This makes this pair of interact ing galaxies 
an ideal target for studying the environments of XRBs. In IClark et al.l (120071 . henceforth 
Paper II) we performed an extensive study of the XRB environments using Chandra X-ray 
images and J and K s IR images (Paper I). Our present paper will expand on our previous 
study by exploring the relationship between XRBs and cluster mass in the Antennae. 

Recent theoretical models of young, massive cluster evolution provide a framework for 
comp arison to our observational study. Two in particular, (|Oskinoval2o"oir 



Sepinskv. Kalogera. fc Belczynsl* 



20051 ). incorporate the fraction of XRBs per cluster in their models. lOskinoval (120051 ) use a 
popu lation synthesis code to study the evolu tion of X-ray emission in young, massive clus- 



ters. ISepinsky. Kalogera. fc Belczynskil (120051 ) investigate the role of supernova kicks in XRB 
expulsion from the parent cluster using the population synthesis code, Star Track. They also 
incorporate the number of XRBs for a range in cluster mass. We will compare our mea- 
surements of the fraction of XRBs per cluster in the Antennae to those predicted by these 
models. 

We organize our paper as follows: In §2 we give a brief summary of our previous work 
on the Antennae and then define a quantity, 77, relating the XRB fraction to cluster mass in 
the Antennae. In our analysis, we estimate cluster mass using K s luminosity, which depends 
non-trivially on the assumed cluster age. We explore cluster age/luminosity relations and 
their impact on our mass estimates in §3. We compare r\ to the measured value predicted 
by theoretical cluster evolutionary models and present conclusions in §4. 
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2. Observations and Data Analysis 

2.1. Infrared Images 

This paper is based on infrared (IR) J (1.25 /im) and K s (2.15 /im) images of the 
Antennae galaxies. We initially presented these data and discussed the details of their 
reduction in Paper I. In summary, 20-minute total exposures in each filter were acquired 



using the Wide-field InfraRed Camera (WIRC - see I Wilson et al.l ( 120031 ) for details) on the 
Palomar 5-m telescope during the night of March 22nd, 2002. In Paper II, we made a frame- 
tie between the IR and X-ray images using IR counterparts to circumvent the poor absolute 
astrometric accuracy of Chandra (~ 1'.'5). We matched seven IR sources from the WIRC 
images with Chandra X-ray point sources. Using a least squares fit of a linear matching 
function we tied Chandra right ascension and declination to WIRC x, y pixel positions. The 
rms positional uncertainty is ~ 0'.'5. With a strong astrometric frame-tie in place we were able 
to accurately identify IR counterparts to X-ray sources. We found 19 IR counterparts within 
1'.'5 of an X-ray source, 13 of which were within l'.'O of an X-ray source. After estimating the 
IR source density, we predict only two of the "strong" matches (separations < l'.'O) and three 
of the "possible" matches (separations between l'.'O - 1'.'5) are due to chance superpositions 
of unrelated objects. 

In Paper II, we pointed out two important implications for these results. First, that 
there is clearly a significant excess of IR counterparts within l'.'O of the X-ray sources - 13, 
where we expect only two in the null hypothesis of no physical counterparts. Even including 
the "possible" counterparts out to 1'.'5, we have a total of 19 counterparts, where we expect 
only five are chance superpositions. Secondly, this implies that for any given "strong" IR 
counterpart, we have a probability of ~ 85% (11/13 with a la uncertainty of 0.30) that the 
association with an X-ray source is real. Even for the "possible" counterparts, the probability 
of true association is ~50%. Therefore, regardless of the physical separations between the X- 
ray sources and their IR counterpart, we are confident that the majority of these associations 
are real. 



We note that of the 19 X-ray so urces with counte rparts, two are the nuclei (IZezas et al. 



2002al ). one is a background quasar (IClark et al.ll2005l ). and two share the same IR counter 



part. Therefore, in this paper we will only consider the 15 IR counterparts (of the original 
19) that are star clusters in the Antennae. 



Found using confidence levels for small number statistics listed in Tables 1 and 2 of lGehrelsl (jl986h 
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2.2. Photometry 

We performed aperture photometry in both the J and K s bands on all 15 IR cluster 
counterparts plus an additional 204 clusters identified by eye in these IR images of the 
Antennae (see also Paper I and the tables there- in). We defined our aperture as ~ 3a of 
the Gaussian PSF, where J had a full width at half maximum (FWHM) of l'/2 and K s had 
a FWHM of Of! 9. We measured a mean and median sky background flux in two separate 
annuli between ~6 - 10cr of the PSF. To account for the exceptionally crowded field of the 
Antennae, we employed the use of background arcs instead of annuli. Multiplying these four 
measurements by the area of the central aperture and then subtracting these from the flux 
in the central aperture yielded four separate source flux measurements. We defined error 
in sky background, a s k y , as the standard deviation of the four measured source fluxes. We 
also considered Poisson noise, a a du, defined as the total source flux divided by the square 
root of the gain for the WIRC instrument. The gain for WIRC during the observations was 
2e~DN~ 1 ( Wilson et al]|2003r fl Adding a sky and a a d u in quadrature, we computed the total 



error in flux, <Jfi ux . We converted fluxes to magnitudes using a bright, 2MASS star in the 
field and defined the error in magnitude, a m , as afi ux divided by the mean flux. Typical 
errors in magnitude were ~0.06 mag in both bands, with no error above 1.0 mag. 

To estimate cluster masses we needed to compute K s luminosity (M Ks ). We computed 
Mk s using reddening derived from (J — K s ) colors (Paper II). Assuming all clusters are 
dominated by O and B stars, their intrinsic (J — K s ) colors are ~0.2 mag. Approximating 
this value as mag, this allowed us to estimate Ak s as ~ ( J — K s ) b s /1.33 using the extinction 



law defined in ICardelli. Clavton. fc Mathisl Jl989T l. 



2.3. XRB-to-Cluster Mass Fraction 

We assume for now that cluster mass is proportional to K s luminosity - i.e. that the 
stellar composition of all clusters is the same. We defined a luminosity cutoff for statistical 
purposes as Mk s = —13.2 mag (see Paper II for details). We binned the data by 0.2 mag 
in Mk 3 and then calculated an average flux per bin. Computing the fraction of the total 
number of clusters per average flux of each bin, we took this as the probability of finding a 
cluster with a specific mass. In Figure 1 we compare this probability for both clusters with 
X-ray sources and all clusters in the Antennae. This shows that XRBs are more common in 



2 At the time of the Antennae observations, WIRC was equipped with a Hawaii-1 lKxlK detector and 
this is the gain for it. 
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more massive clusters. 

This result is not surprising - as star cluster mass increases, so does the number of 
massive stars in it. Through stellar evolution, a certain fraction of these stars will die in 
supernova explosions, leaving behind neutron star or black hole remnants. In turn, a fraction 
of these stellar remnants will retain/acquire a mass-donating companion star, becoming a 
detectable XRB. Thus, through sheer numbers of stars in more massive clusters, we expect 
a greater likelihood of finding XRBs in them. This leads us to two interesting questions: 
1) quantitatively, what cluster mass will more likely produce an XRB and 2) is there some 
intrinsic property of massive cluster physics that favors the production of XRBs, beyond 
simple scaling with mass? 

We believe that our large sample of IR-to-X-ray associations provides the first dataset 
sufficient to estimate the answers to these questions for the Antennae galaxies. In our 
approach to answer these questions we explore the relationship between the number of X- 
ray detections per unit mass as a function of cluster mass in the Antennae. We can formalize 
this expression in the following equation: 



Here, N X (M C ) is the number of detected X-ray sources with an IR cluster counterpart, 
N C i(M c ) is the number of detected clusters, and r](M c ) is the fraction of X-ray sources per 
unit mass, all as a function of cluster mass, M c . 

If rj(M c ) increases or decreases over a range in M c , this means there could be something 
peculiar about massive cluster physics to favor or suppress XRB formation. In contrast, a 
constant i](M c ) across all M c would indicate that more massive clusters are more likely to 
have an XRB simply because they have more stars. 

While rj(M c ) is a powerful tool in studying the number of XRBs per cluster, it requires 
that we know the mass of each star cluster. However, extrapolating the masses of the 
Antennae clusters from our photometric data required models that called for estimates of 
ages and metallicities. While we successfully constrained these inputs and determined cluster 
masses (see below and §3), we first sought to compute i](M c ) in terms of a purely observable 
quantity - flux. Calculating i](M c ) as a function of i^ s -band flux, t](Fk s ), allowed us to 
investigate non-model dependent trends in r](M c ). 

We calculated i](Fk s ) for clusters with a Kg-band luminosity brighter than the -13.2 
mag cutoff using bin sizes of Fk s = 4 x 10 6 in counts (DN) (Figure 2). This bin size was 
small enough to show a trend in i](Fk s ), but large enough to contain at least two clusters 



N x (M c ) = N C i(M c )-r)(M c )-M c 
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with X-ray sources, allowing us to assign error bars to each value of T)(Fk s ). The errors 
plotted on the graph are the measurement uncertainty in the mean value of the four t)(Fk s ) 
added in quadrature with the Poisson uncertainty of the mean t){Fk s ) in each bin. Due to the 
small sample size per bin, we comput ed these uncertainties using the small number statistics 
formulae described in iKeepingl (Il962l ). Figure 2 shows that t)(Fk s ) is roughly consistent with 
a constant value of 5.4 x F^ 1 with an uncertainty of = 1.8 x F^ 1 . 



In essence, r](F K J is comparing two different mass distributions, Nci(F Kg ), the mass 
distribution for all clusters in the Antennae and N x (F Ka )/ F Ka , the mass distribution for 
clusters with X-ray sources, normalized by flux. If there is a constant number of X-ray 
sources per unit cluster mass as suggested by Figure 2, then these two mass distributions 
should be the same. We can further corroborate this result by comparing Nx{Fk 3 ) j F^ s 
and Nci(Fk s ) using a two-sided Kolmogorov-Smirnov (K-S) test. The K-S test yielded 
a D-statistic of 0.75 and a probability of 0.107 that they are related. Considering the 
separate cluster mass populations as two probability distributions, each can be expressed 
as a cumulative distribution. The D-statistic is then the absolute value of the maximum 
difference between each cumulative distribution. This test quantitatively demonstrates that 
there is nothing peculiar about massive clusters in the Antennae with associated XRBs. 
We also computed the Pearson r linear correlation coefficient between Nx{Fk s )/Fk s and 
Nci(FK a ), finding a value of 0.99. Since a value of 1 means a perfect linear fit, this value of 
r further substantiates the observed relationship in t}(Fk„)- 

We then converted t}{Fk s ) into the more conventional units of solar mass. Here we 
assume all clusters a re coeval. Selecting Mk s listed in the Bruzual-Charlot (BC) cluster 
evolutionary mod els (IBruzual fc Chariot 120031 ) for a 20 Myr, 1M cluster as a typical value 



in the Antennae ( Whitmore et al 



19991 ) . we converted the model M Ks to F Ks using the 



standard relationship between luminosity and flux. Multiplying t](Fk s ) by Fk s we converted 
t]{Fk s ) to solar masses: assuming a cluster metallicity of z = 0.02, r\ = 5.8 x 10~ 8 Mq 1 
with an uncertainty of = 1.9 x 10 _8 Mq , while assuming a metallicity of z = 0.05, 
r] = 8.9 x IO^Mq 1 with an uncertainty of = 3.0 x lO^M" 1 . 

While we assume all clusters are ~20 M yr old, we note that the actual range in ages 
should be ~ 1-100 Myr (IWhitmore et al.lll999l ). The BC models indicate that clusters in this 
age range could vary by a factor of as much as 100 in mass for a given K s luminosity. Since 
7] is a function of mass, incorrectly assigning cluster ages has the potential to significantly 
impact rj. In the next section, we explore how differences in cluster age can affect the value 
of rj. 
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3. Effects of Age and Slope in 7/ 

We investigated the effect differences in cluster age has on 77 by assuming three indi- 
vidual age distributions for the Antennae clusters: an instant burst in which all clusters 
are the same age, a uniform dist ribution, and a distribution of the form, dN/dr oc r _1 
( Fall. Chandar. &: Whitmorel 2005). In each case, we assumed all clusters have solar metal- 



licity (z = 0.02) flWhitmore et al. 



1999 ). 



In addition, we address the issue of whether the functional form of f](M c ) can be fit by 
a slope or if, indeed, r)(M c ) is consistent with a single value, by performing a x 2 test. We 
computed the between the 77 values for each bin and the mean value of 77, and computed 
the E% 2 between the 77 values for each bin and a fitted line to these values. The difference 
between the two value of Ex 2 , A£x 2 , indicated the significance of the fitted slope. If AS% 2 
was less than one, with in one S^ 2 deviation, then this indicated that the fitted slope was 
insignificant. While a value of AS% 2 between one and two, with in two S% 2 deviations, 
indicated only a weak slope. Any value of AS% 2 greater than two meant a significant, 
nonzero slope in 77. For our initial case, where we estimated cluster mass using Fk s we found 
AS% 2 = 1.9, suggesting a weak, but non-negligible slope, in 77. 

Considering the case of an instant burst in star formation, we assigned the same age 
to describe all clusters and picked several such values in the range 1 - 100 Myr. Applying 
the BC models, we converted r/(i 7 V s ) to units of solar mass for a range in ages. For each 
distribution in 77 we computed a mean value (Figure 3). The mean for these values is r)i nsta nt 
= 3.3 xl0~ 8 Mq 1 with a standard deviation of a instant = 2.9 xlO -8 Mq 1 . 

Performing our x 2 test for each distribution in 77 (see Figure 3), we found a mean in 
AS% 2 of 0.86, with a standard deviation, ctasx 2 = 1-7- This seems to indicated conflicting 
results. For some cases in assumed cluster age there is no significant slope in 77, while other 
cases show a pronounced slope. This becomes more evident by examining Figure 3, which 
shows a range in AS% 2 of ~0-6.5. Incidentally, for 7/20, AE% 2 = 2.2, while for 77 10 o, AS% 2 
= 4.8xl0~ 3 . 

Next, we assumed a uniform age distribution for our Antennae cluster sample. Picking 
ages at random from a uniform distribution between 1 - 100 Myr, we assigned an age to 
each cluster in our sample. Applying the BC models, we computed each cluster's mass 
based on the assigned age, produced a mass distribution, and then calculated a mean 77. 
Performing a Monte Carlo (MC) simulation, we recreated cluster mass distributions 10,000 
times, producing a large sample of 77's with a mean of Tj uni f orm = 5.4 x 10~ 9 Mq 1 and &rj uniform 
= 6.3 xl(T 10 Mq 1 (Figure 4). 

We then administered our x 2 test for each realized distribution of 77, finding a mean 
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AS^ 2 = 5.9 and ctas x 2 = 5.2. In some cases, rj didn't exhibit a positive slope. Therefore, 
we investigated whether the slope in 77 tended to be mostly positive or mostly negative by 
multiplying each AX% 2 statistic by the sign of the slope. In doing so, we found a mean 
ATjX 2 = 5.8 and 2 = 5.3. Thus almost all distributions in 77 for the uniform case had a 
positive slope. 



A more re alistic approach is assuming the clu ster ages are defined by a power law (PL): 
dN/dr oc r _1 ( jFall. Chandar. fc Whitmordl2005l ). These authors derived their relationship 
using HST UBVIH a observations of ~1 1,000 clusters. Fitting BC models to photometry of 
each cluster, they generated an age distribution. In our analysis we picked ages at random 
according to this distribution. Mirroring the procedure used for the uniform distribution 
case above, we created a sample of 10,000 r] values. We found a mean for this sample of fj PL 
= 8.7xKT 9 M' 1 and a a-, PL = 1.2xl0~ 9 M^ 1 . 

Performing our x 2 test as we did with the uniform age distribution case, we found a 
mean AS% 2 = 5.1 and o"as x 2 = 6.0. Accounting for variations in the sign of the slope in 
77, we found AS% 2 = 4.3 and 2 = 6.6. Therefore, the trend in the slope of 77 remains 



positive, but not as significantly as the uniform age distribution case. 

In a final scenario, we used the cluster ages listed in the electronic table available through 



Mengell (120051 ) to fit ages to 144 c lusters in our sample, including all 15 clusters associated 
with X-ray sources. iMengell (120051 ) derived ages by using three age indicators - UBVI and 
K s broadband photometry to break the age/reddening degeneracy, Ha and Br7 emission to 
identify clusters less than 7 Myr, and CO band-head absorption from narrow-band images 
for clusters ~10 Myr. They then fit these data t o theoretical s pectra for ages < 500 Myr 
using a \ 2 minimization technique (for details, see IMengell 120051 ) . 



Following the met hod discussed in §2.3, we used the BC models to convert M Ks to mass 
for those clusters with IMengell (120051 ) age estimates. Using cluster bins of l.OxlO 7 M in 
mass, we computed four values for 77 (see Figure 5). The errors plotted on the graph are 
uncertainties in the mean value of 77 added in quadrature with the Poisson un certainty in 
each bin. Again, we used the small number statistic formulae in iKeepingl (Il962l ) to compute 
these errors. We found a mean in rj of 2.2x10~ 8 Mq 1 and = 1.2 x lO^M" 1 . Applying 



our x 2 test, we found A£% 2 = 1.0x10 2 , implying no significant slope in 77. 

We summarize our age analysis in Figure 6. Comparing the distributions in 77 for the 
fitted ages, instant burst, uniform and power law age distributions, all values for 77 are within 
a factor of four. Assuming an instant burst of 20 Myr, rj differs by a factor of ten (see Figure 
6) . Excluding specific ages for instant bursts of cluster formation, there is little variation in 
the value of 77. 
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Comparing our \ 2 test f° r each of the four different age assumptions indicates inconsis- 
tent results (see Table 1). In the instant burst, uniform and power law cases, the slope in rj 
varies fro m insignif i cant t o distinctly positive, while rj has no significant slope when the ages 
derived by iMengell (120051 ) are fit to our clusters. Therefore, we can not ignore that rj might 
have a non-zero slope and we will discuss the implications of this in the following section. 



4. Summary and Conclusions 



Through the quantity r], our investigation revealed conflicting results with respect to 
the relationship between observed number of XRBs and cluster mass. We performed a x 2 
test on the slope in 77 for a variety of assumed star cluster age distributions and found some 
cases where the slope was insignificant, while other cases showed a distinctly positive slope. 
No slope indicates the observed number of XRBs per unit mass is independent of cluster 
mass, while a positive slope in 77 suggests more XRBs per unit mass are produced in more 
massive clusters. In the following discussion, we will consider both a constant value in 77 as 
well as a slope in 77, and the implications of each result. 

Initially, we estimated cluster mass by fitting BC spectrophotometric models to cluster 
Mx s , assuming all clusters are ~20 Myr. Recognizing that this method depends on cluster 
age, we explored how different assumptions of a cluster age distribution for the Antennae 
affect 77 and showed that 77 varies by a factor of roughly four; although including individual 
ages for the instant bursts case will increase the variations in 77 to a factor of 10. 

We now proceed by comparing the mean value of 77 for the four different assumed age 

distributions, 77 = 1.7xl0~ 8 Mq 1 , to that predic ted by models of you ng, massive clusters. We 

will co mpare rj to theoretical models discussed in lOskinova J2005h and lSepinsky. Kalogera. fc Belczynskil 
J2005h . 



In a recent study presented in lOskinoval (120051 ). the author modeled X-ray emission from 
young, massive star clusters, assuming a closed system with constant mass, no dynamics and 
all stars are coeval, with cluster metallicities of either z = 0.02 or z = 0.008. These models 
predict ~2-5% of all OB stars in a cluster should produce high mass X-ray binaries (HMXBs) . 
In the models in lOskinoval (120051 ). all clusters are assumed to have masses of M c ; = 10 6 M & 
with stellar masses ranging from 1 - 100 M & . Considering the Salpeter initial mass function 
(IMF) of the form £(M) = M M" 2 - 35 and defining stars with masses > 8 M as "OB stars", 
for our purposes here, we estimated 6% of all stars in the model clusters are OB stars. 
Therefore, 1 - 3xl0~ 3 of all stars in a cluster with an initial mass of 3xl0 6 M (set by 
the Salpeter IMF) should produce an XRB. Since the Salpeter IMF implies there are 7xl0 5 
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stars in a cluster, then these stars should produce 7 - 22xl0 2 XRBs - orders of magnitude 
greater than the ~49 obser ved in the Ante nnae. Expressing 77 as a fraction of XRBs-to- 



cluster mass, the models in lOskinoval (120051 ) suggest rj ranges from 3-7xl0~ 4 M" 1 . These 



values are greater by at least a factor of 1000 from our estimates for r\. Clearly, this predicts 
a much larger n umber of compact object binaries than what we observed in the Antennae. 
Oskinoval (120051 ) note that they were unable to detect HMXBs in three massive (~ 10 4 M & ) 
clusters which they predict should contain between 1-3 HMXBs. If our measured value for 77 
accura tely describes the number of HMXBs formed, then it is not surprising that lOskinova 
( 120051 ) fail to find any. As pointed out by these authors, future modeling of HMXB formation 
is needed to understand the discrepancy between the predictions and observations of XRB 
populations in starburst galaxies. 



In another study, ISepinsky. Kalogera. fc Belczvnskil (120051) use the binary evolution and 
population synthesis code, Star Track ( Belczynski et al. 20021 ) . to investigate the rate of XRB 
formation and ejection from young, massive clusters. This program tracks stellar parameters 
such as radius, luminosity, mass and core mass. The simulations are stopped at the formation 
of a co mpact object. The models include mas s transfer in binaries and include transient 
XRBs. ISepinsky. Kalogera. fc Belczvnskil (120051 ) consider cluster masses ranging from 5 x 10 4 
M & to 5 x 10 6 M & and cluster ages from 1 to ~20 Myr and compute the average number of 
XRBs within 1-1000 pc of the cluster center. 



Considering the typical cluster age in the Antennae is 20 Myr, ISepinsky. Kalogera. fc Belczynski 

( 120051 ) predict a 5 x 10 4 M & cluster should contain 0.13 XRBs, while 15 XRBs should re- 
side in a 5 x 10 6 M & cluster. Here we assume that an XRB is associated with a cluster if 
it is within 100 pc. This separation is equivalent to l'/O at the distance of the Antennae 
(for H =75 km s~ x Mpc~ x ), which is our criteria for an XRB-cluster association (Paper II). 
These model predictions for XRB detections assume a limiting X-ray luminosity of Lx = 
ergs s _1 , but the observed limiting luminosity in the Antennae is 2xl0 37 ergs s~ ] 



5xl0 35 

Using the X-ray luminosity function ( XLF) for the Antennae defined in IZezas fc Fabbiano 



( 120021 ). we scaled the XRB results of ISepinsky. Kalogera. fc Belczvnskil ( 120051 ) to estimate 
what these models would predict for the o bserved number of XRBs in the Antennae clusters. 
Using a XLF power law slope of a = -0.45 (IZezas k, Fabbiano 1120021 ) . the models predict 0.02 



XRBs are observed in a 5 x 10 4 M cluster, while 2.7 XRBs should be seen in a 5 x 10 6 M ( 



© 



cluster, at the luminosity limits of the X-ray observations. Expressing these model results 
as a fraction of XRBs-to-cluster mass, we can directly compare them to ou r measured value 
for 77 in the Antennae. Doing so, ISepinsky. Kalogera. fc Belczvnskil (120051 ) predict 77 ranges 



from 4 — 5 x 10 7 M^ 1 , with in a factor of five from our predictions for 77. As mentioned by 



Sepinsky. Kalogera. &: Belczvnskil (120051 ). several caveats exist for their models including: 1) 
assumed binary fraction of unity which could lead to over estimates of the mean number of 
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XRBs per cluster, 2) the stellar, power-law IMF can affect the XRB fraction per cluster, and 
3) changes in the half-mass radius can strongly influence the median XRB distance from the 
cluster. These factors could potentially explain the discrepancies between their models and 
our observations. 

Since some forms of rj exhibit a distinctly positive slope, these models may not always 
apply to 77. More importantly, a positive slope has implications for star formation scenarios 
in clusters. As mentioned above, a positive slope implies more XRBs per unit mass are 
produced in more massive clusters. If an abnormally large number of XRBs exist in a star 
cluster, then the progenitors of their compact objects should also posses an over abundance. 
Since the progenitors are massive stars, this implies star formati on in massive cluste rs favors 



stars at the heavier end in mass. This is not unusual. Work by IStolte et al.l (120051 ) suggest 



some of the largest clusters in the Milky Way could have a top heavy mass function. 

In this paper we introduced the quantity, 77, relating the fraction of X-ray sources per 
unit mass as a function of cluster mass. Applying this function to the Antennae, we re- 
vealed several important environmental implications for the X-ray sources in the Antennae. 
Specifically, 77 predicts a far differ e nt rel ationship between XRB formation and cluster mass 
than that predicted by "Oskinova) ( 2005 ) and is broadly consistent with that predicted by 



Sepinsky. Kalogera. &: Belczynskil (120051 ). Clearly, future cluster modeling with particular 
emphasis on the relationship between the number of XRBs in a galaxy and the galactic 
cluster environment is essential to explain our current observations. Furthermore, a \ 2 test 
demonstrated the functional form of rj did not always remain consistent with a single value, 
but for some assumptions for a cluster age distributions, the slope had a significantly positive 
value. While this could imply a top heavy mass function in massive clusters, our statistics 
are small. We plan to enlarge our statistical base by extending our observational study to 
additional starburst galaxies. We can then address whether the properties of r/ depend on 
an individual galaxy or are consistent across all galactic environments. 
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Table 1. AS^ 2 Statistics 



Age Test 


A£ x 2 




Median AS^ 2 


flux 


1.9 






20 Myr 


2.2 






100 Myr 


4.8xl0~ 3 






Mengle 


l.OxlO- 2 






Instant 


0.86 


1.7 


0.10 


Uniform 


5.9 


5.2 


4.7 


Uniform 1 


5.8 


5.3 


4.7 


Power Law 


5.1 


6.0 


2.8 


Power Law 1 


4.3 


6.6 


2.1 



Note. - For instant, uniform and power law 
cases, AE% 2 statistics are mean values. 

Multiplied ASx 2 by the sign of the slope for 
the fitted line. 
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Fig. 1. — Here we plot Mk s versus the number of clusters in each bin divided by mean flux 
in that bin. Each bin is 0.2 mag. Arguing that mass is proportional to flux, this graph 
shows the probability of finding a cluster with a given mass. The dashed line signifies the 
magnitude cutoff, M Ks = —13.3 mag. 
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Fig. 2. — This figure displays rj(F Ks ) plotted versus M Kg . The bins are F Kg = 4 x 10 6 DN" 1 
in size. Error bars are the uncertainties in the mean value of rj added in quadrature with 
the Poisson uncertainty in each bin. The dotted line is the mean of the four rj(F Ks ) values. 
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Fig. 3. — Assuming an instant burst of star formation in the Antennae, we plot the mean 
value of r] for a range in ages between 1 - 100 Myr. Notice the factor of ~10 range in rj as 
well as the degeneracy in 77 in this age range. Error bars are uncertainties in 77. 
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Fig. 4. — Comparison between uniform and PL Monte Carlo simulations of rj. The peaks 
of each distribution vary by a factor of ~2 in t], indicating r) does not significantly change 
when we assume different age distributions for the Antennae. We also plot the values of rj 
for instantaneous bursts at three different ages. 
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Fig. 5. — Here rj is plotted versus c luster mass in units of M . In this case, we computed 
cluster mass using ages provided by iMengell (J2005|) (see §3). The bins are M Q = lxlO 7 M 
in size. Error bars are the uncertainties in the mean value of rj added in quadrature with 
the Poisson uncertainty in each bin. The dotted line is the mean value of the four t)(Mq). 
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Fig. 6. — Here we summarize how age affects rj, assuming four different age distributions 
for the An t ennae clusters: instant burst (a), uniform (b), power law (c) and derived ages 
by iMengell (120051 ) (d). Each value is the mean rj and includes 1-a error bars. See text for 
details. Also included is r] for an instant starburst of 20 Myr (e) and 100 Myr (f). 



